Non-universal behavior for aperiodic interactions within a mean-field approximation 
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We study the spin-1/2 Ising model on a Bethe lattice in the mean-field limit, with the interac- 
tion constants following two deterministic aperiodic sequences: Fibonacci or period-doubling ones. 
New algorithms of sequence generation were implemented, which were fundamental in obtaining 
long sequences and, therefore, precise results. We calculate the exact critical temperature for both 
sequences, as well as the critical exponent /3, 7 and S. For the Fibonacci sequence, the exponents 
are classical, while for the period-doubling one they depend on the ratio between the two exchange 
constants. The usual relations between critical exponents are satisfied, within error bars, for the 
period-doubling sequence. Therefore, we show that mean-field-like procedures may lead to nonclas- 
sical critical exponents. 
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Considerable attention has been devoted to the inves- 
tigation of systems displaying inhomogeneous or disor- 
dered interactions [l|. These systems have experimental 
relevance, since many of the materials found in nature 
come with impurities. Moreover, modern techniques are 
able to build materials with controlled composition, such 
that two or more different atoms are combined in a given 
order @. Theoretically, one possible question is to which 
extent the introduction of quenched disorder may affect 
the critical behavior of an originally pure system [l| . The 
answer to this question for random disorder is given by 
the Harris criterion Q : if the pure-system's specific-heat 
exponent, a, is positive (negative), the critical behavior 
of the disordered model is diffferent from (the same as 
for) the pure model. 

The discovery of quasi- crystals Q has motivated the 
introduction of another kind of inhomogeneity in a sys- 
tem: its interactions are now modulated by aperiodic se- 
quences, obtained from, for example, deterministc substi- 
tution rules. These sequences have been a subject of in- 
tense research in recent years [|| : numerical Q as well as 
analytical results [6] have been obtained and the so-called 
Harris-Luck criterion is avalailable [7| . According to this 
criterion, the relevance of a given aperiodic sequence is 
connected to the crossover exponent (f> = l — d a VQ (1 — ui), 
where d a is the dimension in which the aperiodic sequence 
acts, Uq is the correlation lenght's critical exponent of the 
pure model and uj is defined through g ~ N u . Here, g 
is the fluctuation in the number of a given letter of the 
sequence and N is the length of the sequence (see be- 
low). Generally speaking, the greater the fluctuation of 
the sequence the easier to move the critical behavior of 
the system away from the homogeneous case. We can use 
many approaches to address the influence of aperiodicity 
on critical behavior. Mean field is perhaps the most used 



approximation in condensed matter problems, ranging 
from complex fluids Q and superconductivity H to Bose- 
Einstein condensation (lp| and magnetism [llTTl2i]. One 
possible realization of a mean-field picture is the Bethe 
lattice: the critical exponents for homogeneous models 
in this geometry are classical (as for usual mean-field ap- 
proximations), namely: (3 = 1/2, 7 = 1, and 5 = 3, for 
example. In the present case, we arrive at the surprising 
result that aperiodic modulation of exchange constants 
may change the critical exponent of a magnetic system 
within a mean-field approximation. Therefore, the appli- 
cation of this approach seems to be generally useful in 
the study of critical phenomena in complex systems. 

Our goal in this paper is twofold: investigate whether 
there are changes in the universality class of magnetic 
systems with aperiodic interactions in high-dimensional 
systems and determine the existence or not of nonclas- 
sical critical exponents within a mean-field framework. 
To achieve this, we define the Ising model on a Bethe 
lattice, which has proven to be a good approximation to 
the critical behavior in dimension three or above. In this 
geometry, no closed loops are allowed: see Fig. [T] for an 
example of a portion of a Bethe lattice with z = 3. The 
models we study is defined by the Hamiltonians: 



H.l/2 — ~ y JnSiSj, Si — ±1, 



(1) 
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where the sum is over nearest-neighbor spins and J n is 
the exchange constant between sites on generations n and 
n + 1 and may have different values, depending on the 
generations and the corresponding letter in the aperiodic 
sequence (although all interactions between the same pair 
of generations have the same value). 

The exchange constant J n is chosen according to the 
respective letter in the aperiodic sequence. The construc- 
tion of the aperiodic sequence is made in the Bethe lattice 
from the exterior to the interior, i.e, the first letters of 
the sequence correspond to the interactions between gen- 
erations in the exterior of the lattice. The sequences we 
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FIG. 1: Portion of a Bethe lattice with coordination number 
z = 3. Spins are represented by black dots. Generations are 
linked by dashed lines (which do not represent interactions 
between spins) and the interaction constants between different 
generations (continuous lines) may have two different values, 
J a or J B - 



have studied are the following: 

(i) Fibonacci sequence, built from the substitution 
rules: A — » AB, B — ^ A. The first stages of this 
sequence are: A — > AB -> ABA -> ABAAB -> 
ABAABABA. This last finite sequence corre- 
sponds to the following sequence of interaction con- 
stants: J aJbJ aJ aJbJ aJbJ a'-, 

(ii) period-doubling sequence: in this case the substi- 
tution rules are A -> AB and S — ► A A. The 
first letters of the interaction constant sequence are: 
J aJbJ aJ aJ aJbJ aJb- 

The geometrical characteristics of aperiodic sequences 
can be obtained from the substitution matrix Ai, which 
connects the number of letters A and B after one appli- 
cation of the iteration rule, namely: 



N 
N 



(n+l) 
A 

(n+l) 



M 



N 
N 



(n) 
A 

(a) 



(2) 



where is the number of letters £ (A or B) after 

n iterations of the substitution rule. In particular, the 
total number of letters, grows exponentially with 

the number of iterations n: 



~ A?, asn 



(3) 



where Ai is the greater eigenvalue of Ai. The fraction 
of letters A (B), pa (pb) in the infinite word, i.e, after 
n applications of the iteration rules, for large n, is pro- 
portional to the first (second) entry of the eigenvector 
corresponding to the greater eigenvalue of Ai (see Ref. 
IT3I for more details). The fluctuation g^> of a given let- 
ter, say A, is defined as: 
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(4) 



It is possible to show that g ~ N^, where u> — 
In I A2 1/ In Ai (A 2 is the smaller eigenvalue of the substitu- 
tion matrix) . In a linear chain, the fluctuation is greater 
for the period-doubling sequence than for the Fibonacci 
one (see below for a discussion of this point). 

The solution on a Bethe lattice is obtained from a map- 
ping between partial magnetizations of two consecutive 
generations; see [1J] for a detailed calculation of this map- 
ping. In the mean- field limit, z — > 00, J n — > 0, such that 
the product J n = zJ n is finite, and for external mag- 
netic field h, the expression simplifies a great deal and 
we obtain: 



in. 



-i = tanh(_fT n ra n + H), 



(5) 



where m, is the partial magnetization in generation i, 
K n = Jn/ksT, T is the temperature, ks is the Boltz- 
mann constant, and H = h/ksT. Note that we took the 
mean-field limit in order to simplify the equations and 
calculate the critical temperature analitically. 

The critical points of the model studied here are at zero 
external field. The exact critical temperature is given by 
the stability limit of the paramagnetic phase, obtained 
from: 



dm 1 



(6) 



{m}=0 



where m n is the magnetization of spins in generation 
n(n — > 00), mi is the magnetization of spins in genera- 
tion 1 and {m} = means that the derivatives are taken 
in the point where the magnetization of all generations 
from 1 to n — 1 are zero. We then obtain: 



Ja 



(1 



(7) 



where T c is the critical temperature, (1 + r) = Jb/Ja, 
and pb is the fraction of interactions Jb in the ther- 
modynamic limit. As discussed above, this fraction is 
proportional to the second entry of the eigenvector cor- 
responding to the largest eigenvalue of the matrix M. 
[l3| . It is important to know T c exactly (or with a very 
good precision) in order to obtain reliable values for the 
exponent (3. 

The magnetization is not uniform; rather, it follows 
an aperiodic sequence. Therefore, we will calculate and 
discuss the behavior of the mean magnetization, defined 
by: 



_ 1 N 
m = — / 



(8) 



where m-i is the partial magnetization in generation i, 
calculated through Eq. ([5]). This is also the physical 
quantity which would be acessible to experiments. In or- 
der to calculate 7, we define the zero-field susceptibility 
X = drn/dh, related to the mean magnetization. 
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At this point, it is worth recalling that for sys- 
tems such that the length-scaling parameter, in the 
renormalization-group context, cannot assume any value, 
thermodynamic quantities may follow log-periodic scal- 
ing laws [HI, [H, Tri| . The mean magnetization at zero 
external magnetic field, e.g., is written as: 



^Ppog(t)] 



(9) 



where t = (T c — T)/T c is the reduced temperature and 
P[x] is a periodic function of x. Therefore, the logarith- 
mic derivative of m at h = is given by: 



d log Tn 
dlogt 



/3 + P[log(i)], 



(10) 



where P[x] is also a periodic function of x. In an analo- 
gous way: 



and 



d log ui 
dlogh 



dlogx 
dlogt 



l/6 + Q[log(h)}, T = T C , (11) 



= j + R[log(t)], h = 0, (12) 



where Q and R are periodic functions of their arguments. 
Our strategy is to calculate the derivatives, do a best 
fit to a log-periodic function plus a constant and ob- 
tain the value of the respective exponents from this fit. 
Some points are worth stressing: in order to calculate 
this derivative numerically, we need two values of m in 
the thermodynamic limit. This is done in the following 
way. We start with an arbitrary value of mj and iterate 
Eq. ([5]) many times, keeping track of the mean magnetiza- 
tion along the process (the number of iterations we used 
to obtain these results are about 10 8 for the Fibonacci 
sequence and ~ 10 11 for the period-doubling sequence). 
After a transient and a sufficient number of iterations, 
the mean magnetization fluctuates around a monotonic 
trend (see Fig. [2j, which does not depend on the initial 
value mi. From these values, we can extrapolate to the 
thermodynamic limit. This procedure can be done for 
zero or non-zero magnetic field; therefore, we can calcu- 
late m as function of t or h and the susceptibility \- 

In Fig. [2] we depict m n /N for the Fibonacci sequence 
for h = 0, where N is the total number of generations of 
the Bethe lattice or, equivalently, the number of letters 
generated for the sequence, together with the extrapola- 
tion procedure. We take two consecutive pairs of values 
of and make a linear extrapolation of each pair to 
the limit 1/N —> oo. The extrapolated value is taken 
as the mean of the two values obtained and its error is 
estimated as half the interval. This is an overestimation 
of the error but we are sure the true extrapolated value 
is within the interval. We have compared our procedure 
with more formal extrapolations like the so-called VBS 
[HI and BST The extrapolated value is the same, 
although the error estimation is somewhat arbitrary for 



FIG. 2: Magnetization as a function of the number of gener- 
ations of the Bethe lattice or, equivalently, of the number of 
letters of the aperiodic sequence, for the Fibonacci sequence. 
We show a linear extrapolation procedure to obtain the ther- 
modynamic limit for the magnetization and its error, for r — 4 
(see text), 102, 334, 155 generations and log 10 (t) = 1CP 2 ' 5 . 



the VBS algorithm and clearly underestimated for the 
BST one. With this procedure, we were able to calculate 
the derivative in Eqs. (JTUJ) , (fTTj) and (112[) in the thermo- 
dynamic limit. 

A crucial point should be mentioned here: only us- 
ing new algorithms were we able to generate very long 
sequences for the period-doubling case. With these al- 
gorithms, we need to store only N letters, in order to 
use an iV 2 -letter sequence. Long sequences allow for a 
reliable and precise extrapolation to the thermodynamic 
limit (see below) . The algorithm, for the period-doubling 
sequence, works as follows. If an N 2 -letter sequence, with 
N = 2™, n integer, is generated and it is written in lines 
of length N, one over the other, it is easy to note that the 
first N —1 colums have the same letter as in the first line 
of the respective column, while the last column repeats 
the first line if n is even or repeat the first line with A 
and B swapped over if n is odd. Therefore, one needs to 
generate and store only the first line, with N letters, to 
be able to work with an N 2 sequence. This algorithm, 
although with more complicated rules, can be used for 
other sequences, which allow for economy in time and 
space. 

Our results can be summarized in the following fig- 
ures, where we depict the derivatives in Eqs. (fTU)) . (fTTj) 
and (fl"2"|) . These equations tell us that the log-derivatives 
are represented by periodic functions of log(i) or log/i, 
with a mean given by the value of the respective expo- 
nent. For the Fibonacci sequence (Fig. [3]), it is clear 
that the aperiodic modulation does not change the expo- 
nent f3 from its classical value, 1/2. Note the behavior 
of the log-derivative when the size of the Bethe lattice 
(or, equivalently, of the aperiodic sequence) is increased: 
for large values of the reduced temperature t, we are out- 
side of the scaling region and Eq. (fTU|) is not expected to 
be obeyed. When the reduced temperature is decreased, 
the correlation length increases and eventually it grows 
bigger than the number of generations of the Bethe lat- 
tice, leading to the finite-size effect seen in the left-hand 
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s-0 r=4; N=433,494,437 
^< r=0; N=433,494,437 
a-^r=4;N=102,334,155 
v-v i=4; N=63,245,986 





FIG. 3: Log-derivative of the magnetization as a function of 
the reduced temperature for the Fibonacci sequence. We show 
the behavior for different values of r and for different number 
of generations of the Bethe lattice. For comparison, the result 
for the uniform model (r = 0) is also depicted. It is clear that 
j3 is the same for the uniform and aperiodic models. 



FIG. 4: Log-derivative of the magnetization as a function of 
the reduced temperature for the period-doubling sequence, 
with r = 1 and r = 6. The log-periodic behavior is clearly 
identified and the values for f3 are different from the classi- 
cal value 1/2. The continuous line and circles (traced line 
and squares) represent r = 1 (r = 6) and the length of the 
sequenc 



size of the curves (they decrease from the classical value 
1/2). When the size of the lattice is increased, this depar- 
ture from the infinite-size behavior takes place at lower 
values of t, as expected. This finite-size effect is also rep- 
resented by the size of the error bars, as depicted in Figs. 
ID [6] and [Jj (we will comment more on this aspect below). 
The same overall behavior is obtained for 7 and 5, which 
assume their classical values, 1 and 3, respectively. 

The scenario is more interesting for the period- 
doubling sequence. Although the sequences we need to 
generate are bigger than in the Fibonacci case, we were 
able to determine precise values for the exponents (3, 7 
and 5, which depend on the interaction ratio r. In Fig. [4] 
we show the log-derivative of to as a function of log(i) for 
r = 1 and r = 6. The scaling region spans three decades 
and the exponent is clearly different from the classical 
value 1/2: a fitting of the data to a log-periodic function 
leads to j3 = 0.5093(4) and (3 = 0.5664(5) for r = 1 and 
r — 6, respectively. The exponent depends on the ratio 
r, as depicted in Fig. [SJ the overall behavior is the same 
as the one found in Ref. H3, but for a different sequence. 
There is a symmetry with respect to r = 0: the exponent 
(3 is the same for 1 + r = I and 1 + r = l/l, i.e, it is 
invariant with respect to the exchange J a <-> Jb- The 
exponent f3 increases for the aperiodic model, when com- 
pared to its uniform value. This is consistent with the 
following picture: the new (aperiodic) fixed points, in the 
renormalization-group sense, must have a greater value 
of v than the one for the uniform model {vq). Therefore, 
assuming v — 1 jy t to be valid, where yt is the usual tem- 
perature scaling exponent in the renormalization-group 
equations, yt is smaller for the aperiodic model. But 
/3 = (d — yh)/yt (yh is the field scaling exponent) and yu 
is seen to vary very little with r. In fact, our results for 
S (see Fig. given by yh/(d — yh), show that, although 
this exponent depends on r for the period-doubling se- 
quence, it varies much less than [3. So, the variation of 
this exponent is linked almost entirely to yt and it in- 
creases for the aperiodic model (in fact, this picture is 
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FIG. 5: Dependence of the exponent /3 on the parameter r. 
The homogeneous system is represented by r = 0. 



also observed for random disorder). In Fig. [Jj we depict 
d(\ogx)/d(\ogt) as a function of logt, for r = 1 and 6: 
again, it is evident the dependence on r. It is clear, from 
this figure, the need of a second harmonic, to adjust the 
curve for r — 6. The increase of the error bars for small 
values of log(f) in Figs. QJ [5] and is a finite-size effect: 
the smaller t, the greater the correlation length, and the 
size of the Bethe lattice (or, equivalently, the size of the 
aperiodic sequence) must be increased, to reach the limit 
where the infinite-volume behavior is obtained. 

We summarize our results in Table HI where we show 
the values of the three exponents for the period-doubling 
sequence and r — 0, 1 and 6. For r = we show the 
exact results on the Bethe lattice: our numerical values 
agree with them within erros bars. The exponents clearly 
depend on r; nevertheless, a usual relation among the 
three exponents, namely (3 — "f/(5 — 1) is satisfied for 
r = 1 and is just satisfied for r = 6 (see lines 3 and 4 in 
Table!]). 

Preliminary results obtained for the spin-1 Ising model 
lead to exactly the same behavior as for its spin-1/2 coun- 
terpart. As the former model is the Blume-Capel (BC) 
one with zero crystal field A and the latter corresponds 
to the BC model with A = —00, we can infer that our re- 
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FIG. 6: Log-derivative of the magnetization as a function of 
the magnetic field, at T — T c , for the period-doubling se- 
quence, with r = 1 and 6. The exponent 1/8 is different from 
the classical value 1/3 and depends on r The continuous line 
and circles (traced line and squares) represent r = 1 (r = 6) 
and the 
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FIG. 7: Log-derivative of the susceptibility as a function of the 
reduced temperature for the period-doubling sequence, with 
r — 1 and 6. The exponent 7 is different from the classical 
value 1 and depends on r The continuous line and circles 
(traced line and squares) represent r = 1 (r — 6) and the 
length of the sequences is 2 32 . 

suits hold for the BC Hamiltonian and for all values of A 
such that the transition is continuous in the uniform case. 
Also, the results we obtain for the Blume-Capel model are 
not expected to depend on which interaction the disorder 
acts on, according to the following renormalization-group 
reasoning. Even if, initially, only the exchange constants 
follow an aperiodic sequence, the crystal field will not be 
uniform in the coarse-grained system, which is equivalent 
to the original one. 

As a final note, we would like to mention that Eq. @ 
assumes that the relation between the number of letters 
A or B on stage n + 1 of the iteration process depends 
linearly on the number of letters A and B on stage n, 
namely: 

N% +1) = aN^+bN^ 

Ni n+1) = cN^+dN^, [ j 



where a, 6, c, and d are the elements ot the matriz M.. 
The former relation does not hold true for the Bethe lat- 
tice, since letters A on different generations of the lattice 
do not lead to the same number of letters A and B (the 

TABLE I: Values for the exponents ft, 7 and 8 for the period- 
doubling sequence, for r = 0, 1 and 6. In the third and fourth 
lines we show evaluations of /3, calculated from 7 and 8 and 
directly from our data, respectively. Note that the results in 
the first column are exact ones: our numerical values agree 
with them, within error bars. Numbers inside parenthesis are 
errors in the last decimal figures. 





r = 


r = 1 


r — 6 


7 


1 


1.0197(2) 


1.1499(4) 


8 


3 


3.0006(2) 


3.0266(9) 


1 

6-1 


1/2 


0.5097(2) 


0.5674(5) 


P 


1/2 


0.5093(4) 


0.5664(5) 



same is true for the substitution of letters B) on the in- 
flation process. Then, as far as we cannot define the 
substitution matrix for aperiodic sequences on Bethe lat- 
tices, we are not able to determine the exponent u> in 
these cases. 

In summary, we have studied the Ising model on a 
Bethe lattice, in the mean-field limit, such that the ex- 
change constants follow two different aperiodic sequences: 
the Fibonacci and period-doubling ones. The results are 
as follows. For the Fibonacci sequence, the exponent 
are classical, namely /3=1/2,7 = 1 and 5 = 3. For 
the period-doubling sequence the exponents depend on 
the ratio of the interaction constants, but satisfying a 
usual relation among critical exponents for any value of r. 
Therefore, we show that mean-field-like procedures may 
lead to exponents' values different from the classical ones. 
In this particular case, although thermal fluctuations are 
supressed by the mean-field character of the procedure 
we use, geometrical fluctuations present in the sequences 
are treated exactly. Preliminary results for lattices with 
finite coordination number z show the same qualitative 
behavior as for the mean-field limit (z — > 00). 
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